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Kinetics of CO2 Hydrogenation on a K-Promoted Fe Catalyst 

Thomas Riedel and Georg Schaub* 

Engler-Bunte-Institut, Uniuersitat Karlsruhe (TH), Kaiserstrasse 12, 76128 Karlsruhe, Germany 
Ki-Won Jun and Kyu-Wan Lee 

Chemical Technology Division, Korea Research Institute of Chemical Technology, P.O. Box 107, 
Yusong, Taejon 305-600, Korea 

CO2 hydrogenation on a K-promoted Fe catalyst was studied in a fixed-bed microreactor between 
300 and 400 °C, at 1 MPa, and with modified residence times in the range of 0.042-21.4 g-s/ 
cm 3 . For temperatures below 360 °C, organic products almost identical with those found in the 
traditional Fischer-Tropsch reaction with H2/CO were found (paraffins and a-olefins). At 400 
°C, formation of carbon deposited on the catalyst became a major reaction. Concerning the 
mechanism of hydrocarbon formation, the effect of residence time resulted in catalyst particle 
selectivity values for hydrocarbons always higher than zero. This indicates that, besides the 
two-step reaction mechanism via CO, a direct hydrocarbon formation from CO2 can occur in 
principle. With a reaction scheme proposed from these experimental results, a kinetic model 
was developed using integration and regression features of ASPEN PLUS. Calculated values 
for CO2 conversion and CO and total hydrocarbon selectivities agree with the experimental data 
within a range of error less than 15%. 



1. Introduction 

Increasing global energy demands during the past 
decades has caused increases in man-made CO2 emis- 
sions. Although correlation with global temperature- 
effects is not clear, efforts are taken to reduce CO2 
emissions, as was claimed, e.g., at the Third Conference 
of the Parties (COP-3) in Kyoto in 1997. 1 Besides 
improving the efficiencies of energy conversion (e.g., in 
power plants) or energy utilization processes, strategies 
include secondary measures, e.g., capture, storage, and 
fixation of CO2 produced. 2 In this context, chemical 
utilization of CO2 in the conversion of H2 from CO2 
neutral sources to liquid fuels via Fischer-Tropsch (FT) 
synthesis is a potential way. 3 

Catalytic hydrocarbon synthesis from synthesis gas 
CO/H2 recently gained special attention under the 
catchword "Gas-to-Liquids" as a way for producing 
liquid fuels from natural gas. 4 Although there has been 
some progress in understanding the chemistry of CO2 
hydrogenation and in identifying potential catalysts, 
information about the kinetics of this reaction is not 
available. This information could be used for feasibility 
studies regarding the CO2 utilization potential of this 
method. 

Hydrogenation of CO2 so far has been studied mainly 
on catalysts used traditionally for the FT reactions. On 
silica-supported Ni and Ru catalysts, the main product 
is methane; only minor amounts of higher hydrocarbons 
were observed under the reaction conditions applied. 5 
For cobalt as an active metal, Riedel et al. 6 studied the 
effect of varying CO/CO2 ratios in the synthesis gas. 
With increasing CO2 content, the product composition 
shifted from a typical FT-type (mainly hydrocarbons) 
to almost exclusively methane. Generally, Co catalysts 
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exhibit no significant activity for the CO shift reaction 
at moderate reaction temperatures 7 or for the formation 
of CO via the C0 2 shift reaction. Therefore, the FT 
regime of specific inhibition of product desorption 
through strong reversible CO adsorption as a prereq- 
uisite on cobalt for the formation of long-chain hydro- 
carbons cannot be obtained under such conditions. 8 

So far, iron exhibits the most promising properties 
for hydrogenation of CO2 to form long-chain hydrocar- 
bons. 9 It is well-known as an active metal in commercial 
FT synthesis (e.g., in the ARGE catalyst 10 ) and is also 
widely used in commercial CO shift reactors, 11 an aspect 
which will be addressed below. As a consequence, much 
attention has been paid in the past on improving the 
catalytic properties of Fe-based catalysts through the 
use of different support materials 12 or the addition of 
promoters, e.g., Mn, Zn. 1314 The most promising Fe 
catalysts for CO2 hydrogenation were obtained using 
AI2O3 as the structural support and potassium as the 
chemical promoter, applied in concentrations of up to 
0.5 mol-K/mol of Fe. 15 - 16 

Previous studies indicate that hydrogenation of CO2 
proceeds via a two-step reaction mechanism. 5,1317 In the 
first step, CO2 is converted to CO in the CO2 shift 
reaction (reverse water gas shift reaction, eq 1). The CO 

C0 2 + H 2 — CO + H 2 0 A R tf° 300 . c = 38 kJ/mol 

(1) 

produced in this way reacts subsequently in the FT 
reaction to form organic products (eq 2). Using a H2/ 

CO + 2H 2 - (CH 2 ) + H 2 0 

a r#°3oo-c = _166 kJ/mol (2) 

CO and a H2/CO2 synthesis gas, respectively, Schulz et 
al. 18 compared the changes in catalyst activity and 
product selectivity during the initial formation/activa- 
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Table 1. Published Values for Apparent Activation 
Energies of the CO2 Shift Reaction (with No Information 
Given about Rate Constant and Inhibition Parameters) 



catalyst 



7Y°C 



EA/(kJ/mol) 



ref 



FeSi 
FeSi-H 
FeSi-K 
Fe 3 0< 



250-400 
250-400 
250-400 
327-427 



Fe 3 04/Cr 2 0 3 327-427 



56, 72 
81 

68, 98 
-80 
-80 



Kaspar et al. 36 
Kaspar et al. 36 
- Kaspar et al. 36 
Spencer 20 
Spencer 20 



tion of a K-promoted iron catalyst. The transient 
behavior observed was nearly identical, comprising a 
sequence of kinetic regimes (episodes); however, the 
duration of the different episodes lasted much longer 
with CO2. At steady-state conditions, almost the same 
product composition was obtained with both synthesis 
gases, indicating that formation of organic products 
from CO2 proceeds with CO as the intermediate. 

Recently, Fiato et al 19 proposed by using laser- 
generated iron carbide catalysts that there is also some 
contribution possible from the direct hydrogenation of 
CO2 (eq 3). This reaction is proposed to proceed via 

C0 2 + 3H 2 -(CH 2 ) + 2H 2 0 

Ar^soo-c = -128 kJ/mol (3) 

dissociative adsorption of CO2 followed by hydrogena- 
tion of the adsorbed carbon species. However, the extent 
to which this reaction contributes to the overall forma- 
tion of organic products remains open at the moment. 

Concerning the kinetics of the CO2 shift reaction, only 
limited information is available. Published experimental 
data indicate that the reaction is very fast over several 
metallic and oxidic catalysts, having a rather low 
activation energy (Table 1) and high reaction rates, in 
comparison with the CO shift reaction (the reverse 
reaction). 20 For the CO shift reaction, tshi, a kinetic rate 
equation such as that shown in eq 4 is typically used. 21 



r SHl ~ *SH1 



PC0PH 2 0 ~ PcO^H/^eq 
PCO + a SHl,H 2 0PH 2 0 



(4) 



Here, the effect of chemical equilibrium limitations is 
considered by introducing the equilibrium constant K\ 
(eq 5). Summaries of rate equations for the CO shift 



log^ q = log(^^) 

\PH,0PC0/ eq 



2073 
77K 



- 2.029 (5) 



reaction were published, e.g., by Zimmerman and 
Bukur. 22 

For the CO2 shift reaction, no information is available 
about rate equations. Therefore, an expression analo- 
gous to the one for the CO shift reaction was used as 
the starting point in this study. However, a term 
accounting for a possible inhibition by CO2 was added 
in the denominator (Table 2). 

In the FT reaction (eq 2), Dry 23 found that, for a 
reduced Fe catalyst at 225-265 °C at low conversion 
and low partial pressure of H 2 0, the rate of total 
synthesis gas conversion is proportional to the partial 
pressure of H 2 . At higher conversion, i.e., when the 
partial pressures of the products are significant, a more 
general rate equation should be used as a starting point 
for kinetic modeling, 24 which is based on the one 
proposed by Anderson 25 (eq 6). To this equation was 



Table 2. General Form of Rate Equations for C0 2 Shift 
Reaction (rsi!2), FT Reaction (rpr), and Direct CO2 
Hydrogenation (n>n) Used in the Present Study 



reaction 



rate equation 



CO2 shift r SH2 — & SH2 - 

Fischer- Tropsch = k F 



PcOjPHa PC0PH 2 c/^eq 
l PC0 + a SH2.H 2 oPH a O + bsm.CO./'cO., 

PcoPhj 



direct 

hydrogenation 



T PC0 + a FT.H 2 oPH 2 0 + frFT.CO./cO.j 
PC0./>H 7 



PCO + °DH.H. 2 0Ph 2 0 + ^DH.COrPcO.j 



added in the present study a term for possible inhibiting 
effects of C0 2 (Table 2). 



PCO + a FT.H 2 oPH 2 0 



(6> 



Other rate equations were proposed, for example, by 
van Steen and Schulz 26 and Huff and Satterfield, 27 and 
summarized by Zimmerman and Bukur 22 and Espinoza 
et al. 28 In all of these equation^, a term for the inhibiting 
effect of water vapor is included. Only Ledacowicz 29 and 
Nettelhoff, 30 both from the same research group, con- 
sidered exclusively the inhibition by CO2 instead of H 2 0. 
Zimmerman and Bukur 22 studied FT kinetics of two 
different Fe-based catalysts under identical reaction 
conditions. Considering the inhibition of both H2O and 
CO2 gave the best results in the case of a precipitated 
Fe— Cu-K catalyst, whereas inhibition by CO2 was 
negligible in the case of a commercial Fe catalyst 
(Ruhrchemie LP 33/81). Table 3 gives an overview of 
the kinetic parameters for the FT reaction which can 
be found in these publications. 

Espinoza et al. 28 pointed out that kinetic studies 
should be performed on the specific catalyst chosen. In 
addition, they should also cover the proposed range of 
operating conditions for design and scale-up studies. 
This is clearly shown by the variety of rate equations, 
especially the inhibiting effect of H 2 0 and 002, which 
indicates a strong link between the composition of the 
catalysts, the reaction conditions applied, and the final 
rate equation and kinetic parameters obtained. 10 For 
FT catalysts, highly active for the CO shift reaction (e.g., 
with K as the promoter 11 ), leading to low H 2 0 and high 
CO2 partial pressures, inhibition of C0 2 (and H2O) 
should be considered, whereas at high H2O/CO2 ratios, 
the inhibiting effect of CO2 is minor compared to the 
effect of water. These findings suggest that in the case 
of CO2 hydrogenation an inhibiting effect of H 2 0 as well 
as of C0 2 has to be considered in the rate equations for 
the FT and CO2 shift reactions. Because no information 
is available about the kinetics of the direct formation 
of organic products from CO2, a rate equation similar 
to that of the FT-type appears to be appropriate (Table 
2). 

2. Experimental Section and Data Analysis 

2.1. Experimental Procedures. The Fe catalyst 
used in the present study was prepared through quick 
precipitation from hot aqueous nitrate solutions of the 
metals Fe, Al, and Cu, with an aqueous NH3 solution. 
The washed, dried (12 h at 120 °C in air), and crushed 
catalyst was then impregnated with a K2CO3 solution. 
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«F"Flniol/(s g*Mra)| 


/S/yVKJ/mOl ) 


ret 


225-265 


na 


0 


na 


63 


Dry 23 


270 


4.51 


0 


64.8 


(89)° 


NettelhofTet a!. 30 


240 


0 


0.19 


36 


(sir 


NettelhofTet al. 30 


250 


4.8 


0.33 


288 


(86)° 


Zimmerman and Bukur 22 


250 


4.5 


0 


136.8 


my 


Zimmermann and Bukur 22 


300 


17.5 h 


0 


0.396 


66.5 C 


Liu et al. 37 


270 


3.016 


0 


1432.8* 


na 


Raje and Davis 38 



Hj/CCVAr 




organic 
proOod* 



hot wax (rap 



ampoule sampling 
device 



Table 3. Published Values for the FT Rate Equation As Used in the Present Study (Table 2), Experimentally 
Determined with H*/CO Synthesis Gases 

catalyst 

Sasol Fe catalyst 
precipitated Fe 
fused Fe BASF S6-10 
100/0.3/0.2 Fe/Cu/K 
Fe-based Ruhrchemie LP 3/81 
coprecipitated Fe-Cu-K 
100/0:015/0.005 Fe/Si/K 

0 Apparent activation energy for first-order kinetic rate constant [T = 220-265 °C (Zimmerman and Bukur 22 ) and T = 230-250 °C 
(NettelhofTet al. 30 )!. b Heat of adsorption determined to be 65.4 kj/mol. c Preexponential factor *ft.o determined to be 4.38 x 10 5 mol/ 
(s-g-MPa); T = 220-300 °C. d k related to the mass of active metal Fe. 

heating with 2 K/min to 400 °C, and holding time 10 
h). The reactor was then cooled to the first reaction 
temperature of 300 °C, the reaction pressure of 10 bar < 
was adjusted, and the desired reactant gas flow was set. 
The composition of the gas was H^CCVAr = 71.76/24.1/ 
4.14 (mol %); Ar was used as the internal standard for 
the gas chromatographic (GO analysis. Initial activa- 
tion/formation of the catalyst was monitored during a 
period of several days. Then, after steady-state condi- 
tions were ensured for 24 h (also done for every further 
data point), new reaction conditions were adjusted. 
When the amount of catalyst (0.1-2.5 g) and the gas 
flow rate (7-140 cnrVmin) were changed, the modified 
residence time T m0£ j was varied between 0.042 and 21.4 
g-s/cm 3 . After the reproducibility of the experiment was 
checked, the next temperature, 330, 360, and 400 °C, 
was set. Analysis of the gases H2, CO2, CO, CH 4 , and 
the internal standard Ar was done using an on-line GC 
with a thermal conductivity detector (TCD). Organic 
products were analyzed using an ampule sampling 
technique and adapted GC with a flame ionization 
detector (FID). 31 The carbon balance was closed within 
-2 and + 5 C % (except at 400 °C as discussed in section 
4.1). 

2.2. Data Analysis and Definitions. With the data 
obtained from the GC-TCD analysis, the ratio of the 
peak areas of the inorganic gases and of methane 
related to the peak area of the internal standard Ar was 
calculated. When these ratios were multiplied with 
earlier determined calibration factors and the exact 
known molar flow of Ar, the molar flows of these gases 
were obtained, which were used to calculate the conver- 
sion of CO2 and the yield of CO and CH 4 . 

The peak areas of the organic products obtained from 
the GC-FID analysis were referred to the peak area of 
methane. With the molar flow of CH 4 known from the 
GC-TCD analysis, the molar flows of the organic 
products were calculated and used to determine' their 
yields and selectivities on a carbon basis (in C %). A 
semilogarithmic plot of the mole fraction of organic 
products against the carbon number gives the Anderson- 
Schulz-Flory (ASF) plot, 11 where the slope yields the 
chain growth probability a. 

Kinetic rate parameters of CO2 hydrogenation were 
obtained by varying the modified residence time ex- 
pressed as r m0 d (eq 7). 



Figure 1. Flow scheme of a fixed-bed microreactor (BPR: back 
pressure regulator). 

Finally, the catalyst was dried again and crushed to a 
particle size below 100 ^m to avoid intraparticle mass- 
transfer limitations, even in the catalyst pores filled 
with liquid FT products. The desired composition of 
the catalyst was 100/13/10/10 Fe/Al 2 03/Cu/K (mass 
ratio). Atomic absorption spectroscopy analysis showed 
a final composition of 100/13.3/10.7/8.9 Fe/Al^OyCu/K. 
Brunauer-Emmett-Teller surface area was deter- 
mined as 151 m 2 /g for the dried catalyst, which has a 
bed density of 1.15 g/cm 3 . 

The relatively high concentration of K was chosen 
because alkali promotion appears beneficial for the CO2 
hydrogenation reaction. 12 However, too large amounts 
of K are reported to cause deactivation of the catalyst. 18 
AI2O3 was used as the structural promoter because it 
improves iron dispersion, and copper was added to ease 
the reduction of the catalyst. 

CO2 hydrogenation was studied in a fixed-bed mi- 
croreactor (flow scheme in Figure 1). The reactor tube 
(i.d. = 0.9 cm, / = 32 cm) was filled with the desired 
amounts of catalyst diluted with silica particles {dp = 
250-400 ^m) in a weight ratio of catalyst to silica of 
1/10. The catalyst particles covered the larger silica 
particles as an adhering layer; thus, a uniform gas flow 
at low pressure drop could be achieved, and the high 
degree of catalyst dilution with the inert material 
allowed for operation at almost isothermal conditions. 
The difference in temperature along the reactor middle 
axis measured with a moveable thermocouple was in 
the range between ±1 and ±3 °C, depending on the 
reaction temperature applied and the length of the 
catalyst bed used (8-15 cm). 

After the reactor was filled, the dried catalyst was 
calcined in Ar (40 cm 3 /min, heating with 2 K/min to 400 
°C, and holding time 10 h) and. subsequently reduced 
at atmospheric pressure in a H2/Ar flow (5% H2 in Ar, 
5 cm 3 (0 °C, 0.1 MPa) of H 2 per min and per g of Fe, 



^mod 



V. (0 °C, 0.1 MPa) 



(7) 



Concerning the reactor model to be applied for the 
kinetic analysis of the experimental data, preliminary 
theoretical studies suggest that, under the reaction 
conditions applied, axial dispersion (Pe nx ~ 2, l^d p > 
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400, resulting in Bo ax > 800) and wall effects {djdp > 
10) can be neglected. Because of the small particle size 
applied, internal mass-transfer limitations are not 
significant. 18 Although a limited variation in tempera- 
ture was observed along the reactor axis (max ± 3 °C 
at a reaction temperature of 400 °C), the assumption of 
an isothermal ideal plug-flow reactor (PFR) was made. 
As a result, a steady-state plug-flow model is used for 
each component i (eq 8). F, is the axial molar flux in 



— = Z Wc 

ax i 



(8) 



mol/(cm 2 -s), x, the reactor axial length in cm, v iJt the 
stoichiometric factor of i in reaction./', r,-, the correspond- 
ing reaction rate in mol/(g-s), and p cat , the catalyst bed 
density in g/cm 3 . 

Reaction rates are defined on the basis of the total 
mass of the catalyst used (eq 9). In this equation, £ 
represents the extent of each reaction. 



= J_ 



^1 = 
dt 



l l./dnA 



(9) 



3. Numerical Methods 

Thermodynamic calculations and kinetic analysis of 
the experimental data obtained were done using ASPEN 
PLUS. 32 Equilibrium conversion of C0 2 for a given 
synthesis gas composition was calculated by minimizing 
the Gibbs free energy, ArG, for a given set of compo- 
nents allowed to be present at equilibrium. For this 
purpose, the reactor model RGIBBS in ASPEN PLUS 
was used. 

For kinetic analysis, i.e., determination of kinetic 
parameter values for CO2 hydrogenation, the reactor 
model RPLUG was chosen, which provides a rigorous 
simulation of an ideal PFR. 

Kinetic parameter values were determined via mini- 
mization of the sum of square errors f min (eq 10), with 
7v*ex P , the number of experiments (3, each for one 
temperature), and Ndatn, the number of data points 
obtained for varying r in the each experiment (14, 9, 
and .7, respectively). As variable <l>, the molar flows of 



u-Z I- 



-f 



I 



C |>CXP _ tfecalcl 



'out / 



(10) 



reactor.out/ 



H2 and CO2 at the reactor inlet and of carbon-containing 
components CO2, CO, and C 3 H 8 at the reactor outlet 
were chosen. For the measured data at the reactor inlet 
and outlet, a standard deviation o reflecting the experi- 
mental errors had to be specified. To fit the calculated 
data to the experimental data with the regression model 
DATA FIT from ASPEN PLUS, parameters of the 
reaction rate equations used and of molar flows at the 
reactor inlet were varied within a specified range of 
values, until function /inim calculated with the corre- 
sponding flows at the reactor outlet, was minimized. 32 
Figure 3 shows the result for a reaction temperature of 
300 °C. 




200 



250 



400 



450 



300 350 
T. °C 

Figure 2. Calculated equilibrium conversion of CO2 (curves) 
for an atomic ratio of H/C/O = 6/1/2 |(pn/pco 2 )in = 31 as a function 
of temperature for p = 1 MPa considering the following com- 
ponents: (A) H 2 , C0 2 . CO, and H 2 0; (B) H 2 , C0 2 , CO, H 2 0, 
CH 4 , and C3H8 (representing hydrocarbons); (symbols) experi- 
mental data. Catalyst: 100/13/10/10 Fe/Al 2 On/Cu/K. p = 1 MPa, 
(pH-/pco,)in = 3, and r mo d = 21.4 g-s/cm 3 . 

4. Experimental Results and Discussion 

4.1. CO2 Conversion, CO, and Total Hydrocarbon 
Selectivity. According to the main reaction pathway, 
CO is considered to be first formed via the CO2 shift 
reaction and subsequently converted to hydrocarbons. 
Gibbs free energy calculations indicate that, for the CO2 
shift reaction (eq 1) alone, chemical equilibrium con- 
straints exist in the temperature range of interest here. 
Equilibrium conversion of CO2 increases with increasing 
temperature; however, only about 23% at 300 °C and 
still less than 45% at 450 °C can be converted to CO 
(Figure 2). Considering also the formation of hydrocar- 
bons (CH 4 as the most stable one and C3H8 representing 
a relatively stable species of the organic components), 
a significantly higher equilibrium conversion is possible, 
however, being still far away from total conversion. 
Consequently, a limitation of CO2 shift conversion 
achievable from thermodynamic considerations is obvi- 
ous. This clearly shows the disadvantage of CO2 hydro- 
genation as compared to CO hydrogenation (traditional 
FT reaction), where equilibrium values of Xco are 
practically 100% in the whole temperature range. 

Measured CO2 conversion at the highest value of z mo 6 
applied (ranging from 0.042 to 21.4 g-s/cm 3 ) increased 
from about 45% at 300 °C to almost 60% at 400 °C 
(Figure 2). These values are much higher than the 
calculated equilibrium conversion obtained for the CO2 
shift reaction alone (Figure 2, curve A). 

Measured conversion values of CO2 were found to rise 
sharply with r mo d for low residence times, followed by a 
less pronounced increase at longer residence times 
(Figure 3). Almost doubling of r mot i at high values from 
12 to 21.4 g-s/cm 3 at 300 °C resulted only in an increase 
of CO2 conversion from about 40 to about 45%. This slow 
increase is possibly due to (a) an approach to the 
equilibrium conversion of CO2 (Figure 2, curve B), (b) 
the decreasing partial pressures of the reactants CO2 
and H2, and (c) a possible kinetic inhibition by water 
vapor formed during the reaction. This inhibiting effect 
of water vapor on an Fe-based catalyst is also known 
from traditional FT synthesis starting with H2/CO, 
where it plays an important role in process design and 
scale-up strategies, because product inhibition limits the 
achievable single-pass conversion. 33 Analogous results 
were obtained for the other reaction temperatures 
applied, where with increasing temperature the mea- 
sured CO2 conversion values approached more and more 
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Figure 3. CO2 conversion and CO and total hydrocarbon (HC) 
yields versus modified residence time. Catalyst: 100/13/10/10 Fe/ 
AI2O3/CU/K. T = 300 e C,p = 1 MPa, (p H /pco,)in = 3 (data points, 
experimental; curves, calculated with kinetic model). 



100 








80 

^ 60 
O 




K 




20 
0 


\ • • 


/ Sac 

A • • • 





0 10 20 30 40 50 



Xfj02 . % 

Figure 4. Total hydrocarbon (HC), CO, and "C-loss" (S AC ) 
selectivity versus C0 2 conversion. Catalyst: 100/13/10/10 Fe/Al 2 CV 
Cu/K. T = 300 °C,p = 1 MPa, (pH./pco a )in = 3, r mod = 0.042-21.4 
g-s/cm 3 (data points, experimental; curves, calculated with kinetic 
model; dashed curves, extrapolation of experimental data to 
Xco, - 0). 

the equilibrium conversion. At 400 °C, about 80% of the 
equilibrium conversion was achieved. 

The measured yield values of CO as a function of r mo d 
exhibit the typical behavior of an intermediate product 
in a consecutive reaction (Figure 3). At short residence 
times, a strong increase of Yqo is observed, correspond- 
ing to the CO2 conversion, followed by decreasing values 
at higher residence times when the consecutive reaction 
products, i.e., hydrocarbons, are formed as indicated by 
the curve obtained for the yield of hydrocarbons. The 
same pattern was found for all reaction temperatures. 
These results can be seen as evidence for the proposed 
two-step reaction mechanism for the formation of hy- 
drocarbons from CO2. 513 17 

However, at the lowest CO2 conversion obtained, 
respectively the lowest residence time applied, hydro- 
carbons were also found as products. For Shc, as a 
measure for particle selectivity, extrapolation toward 
zero C0 2 conversion (Figure 4) resulted in a value 
higher than zero, indicating that also hydrocarbons can 
be formed as a primary product, however, in signifi- 
cantly smaller quantities than CO (5.9 C % versus 93.7 
C % at Xco 2 = 3.35%). This result is in agreement with 
Fiato et al., 19 who proposed a possible direct formation 
of hydrocarbons from CO2. Sac in Figure 4, the "C-loss" 
selectivity (corresponding to adsorbed C, carbid, carbon, 
chemisorbed organic species, and nonvolatile organic 
components)., represents a sink term in the carbon 
balance. 

The selectivity data obtained for a reaction temper- 
ature of 400 °C differ significantly from those obtained 
at lower temperatures. This can be explained in terms 
of the M C-loss" selectivity S AC - Below 400 °C, values of 




(paraffins, olefins 
oxygenates) 

SHI ,2 : C0 2 + H 2 CO * H 2 0 

FT : n CO ♦ 2n H 2 *~ C n H 2n + n H 2 0 

DH : n C0 2 ♦ 3n H 2 — »- CnH^ + 2n H 2 0 

Figure 5. Reaction scheme proposed for CO2 hydrogenation on a 
K-promoted Fe catalyst, rsm: reverse C0 2 shift reaction (CO shift 
reaction). rs\\2' CO2 shift reaction, r^j: Fischer-Tropsch reaction. 
tdh: direct formation of hydrocarbons from CO2 (rpr, ^dh: forma- 
tion of olefins). 

Sac were found in the range between -2 and +5 C % 
(e.g., Figure 4 for 300 °C). However, for 400 °C the value 
increased drastically to 15-25 C %. The rate of carbon 
deposition according to the Boudouard reaction (eq 11) 
was found to be proportional to pcdpw^ and to increase 
strongly with increasing temperature. 10 

CO + CO — C + C0 2 A R tf°3 00 . c = -173.5 kJ/mol 

(11) 

At longer residence times, partial pressures of CO at 
400 °C were about 2-3 times that for the lower 
temperatures. Therefore, because no liquid products 
were formed in the reaction, this "C-loss" selectivity 
should be almost exclusively attributable to carbon 
deposited on the catalyst. This is not necessarily in 
conflict with the fact that the catalyst was still active 
at the end of the experiments. Dry 34 reported in this 
context about Fe catalysts used in high-temperature FT 
synthesis that "a common oversimplification found in 
the literature is that carbon deposition necessarily 
resulted in activity loss". After several weeks of opera- 
tion, there can be twice as much carbon as iron in the 
catalyst and the activity is still high. In the extreme, 
there is a "carbonrsupported iron catalyst" in the end, 
with the active iron sites still being readily accessible 
for the reactants. However, it is clear that secondary 
reactions at high temperatures can lead to carbon-rich 
aromatics ("coke"), which results in activity loss. 

Because liquid products and olefins are the desired 
products of CO2 hydrogenation, the high S A c values, 
found for 400 °C, clearly show that temperatures above 
360 °C should be avoided. It also means that for kinetic 
modeling an additional reaction for carbon formation 
must be considered. Therefore, because carbon is an 
undesired byproduct, only the results obtained for 
reaction temperatures of 300, 330, and 360 °C will be 
used to model the formation of hydrocarbons from CO2 
and H2. As a final conclusion from the experimental 
results so far, the reaction scheme shown in Figure 5 is 
proposed for the reaction conditions applied in this study 
and will be used for the kinetic modeling of the reaction. 

4.2. Organic Product Distribution. Among the 
organic products, methane is thermodynamically most 
favored. However, it is also the most undesirable 
product in C0 2 hydrogenation. Therefore, the main 
objective must be to obtain the lowest possible methane 
selectivity. Total S C n 4 values at 300 °C increased with 
increasing residence time parallel to the increasing yield 
of hydrocarbons from 2.5 C % up to a value of about 10 
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Figure 6. ASF plot of volatile organic products (left) and molar olefin content in linear hydrocarbons (right) versus carbon number from 
C0 2 hydrogenation with effects of reaction temperature (X COj = 46-52%, r mod = 15-21.4 g-s/cm 3 ). Catalyst: 100/13/10/10 Fe/AJzOyCu/K 
p = 1 MPa, (p Hj /pco,)in = 3. 



C %. With increasing temperature up to 360 °C, only a 
slight increase of S C h 4 up to 13.5 C % was observed. 
However, at 400 °C, a significant increase was found; 
almost 60 mol % of volatile hydrocarbons were CH 4 
compared to about 30-40 mol % for the lower temper- 
atures at similar C0 2 conversion values of 46-53%. This 
result together with the high tendency for carbon 
formation (S AC ^ 15-25 C %) indicates that at 400 °C 
the FT regime of long-chain hydrocarbon formation does 
not prevail anymore. 

The total organic product distribution can be de- 
scribed in a simplified way by means of the ASF plot. 
At 300 °C, the value of the chain growth probability cii 
(C4-C7) of about 0.56 was found to be almost indepen- 
dent of the residence time and the CO2 conversion. In 
contrast, a slight decrease was found for a 2 (Ci 0 -Ci 4 ), 
with a value of 0.65 for the longest residence time 
applied. Although a small effect of residence time was 
observed, it must be mentioned that the organic prod- 
ucts with N c > 12 amount to less than 3 mol % of the 
organic products found for the conditions applied here. 
With the exception of the methane selectivity, the 
hydrocarbon selectivities and the ASF distribution were 
only slightly affected by variations in the temperature 
in the range from 300 to 360 °C (Figure 6, left). 

From the traditional FT synthesis starting with H2/ 
CO, it is known that linear a-olefins and linear paraffins 
are the main primary products. An adsorbed surface 
species can form either a paraffin in the gas phase via 
associative desorption with adsorbed H or an a-olefin 
via dissociative desorption, leaving an adsorbed H on 
the catalyst surface. 8 The primary olefin selectivity is 
about 80%, but a-olefins can undergo secondary reac- 
tions such as hydrogenation, double bond shift, and re- 
incorporation in chain growth reactions. 35 The experi- 
mentally determined olefin selectivity at 300 °C shows 
a significant decrease for C 2 with increasing residence 
time, reflecting the high reactivity of ethylene for 
secondary reactions and only a moderate decrease for 
the longer hydrocarbons. At the longest residence time 
applied, more than 80% of the products up to C12 were 
olefins. Among these olefins, more than 90% were 
a-olefins. With increasing reaction temperature the 
olefin content increased slightly (Figure 6, right), an 
effect which is known from CO hydrogenation, 35 and still 
more than 90% of the olefins were primary olefins. 
Besides liquid-saturated hydrocarbon products, also 
a-olefins are valuable products, which can be used as 
feedstock for the petrochemical industry. In this respect, 
the catalytic system investigated in this study exhibits 
excellent properties with its high a-olefin selectivity. 



At 300 °C, the selectivity of oxygenates for C3+ was 
found to be in the range of 5-10%. Almost no effects 
were observed for varying residence time or reaction 
temperature. However, within the C 2 fraction, ethanol 
and acetic aldehyde represented about 50% of the 
products formed at 300 °C. Their amount was only 
slightly affected by variations in residence time; how- 
ever, with increasing reaction temperature, the oxygen- 
ates content within the C 2 fraction decreased to about 
30% at 360 °C. 

5. Kinetic Analysis and Discussion 

For the kinetic analysis of the C0 2 hydrogenation 
data, the model of an ideal PFR operating at isothermal 
conditions was applied (see section 2). Because of the 
complex product composition, propane, which has a 
molecular weight similar to the average of the organic 
products formed in the reaction (Nc,av = 2.5-2.8), was 
used as the model hydrocarbon for the stoichiometry of 
the reaction. For the reaction scheme proposed from the 
experimental results (Figure 5), the following set of 
differential equations is obtained: 

C0 2 -j^- = -r S H 2 Pcat - r DH p cat ( 12) 



CO = r SH 2 Pcat ~ ^Pcat ( 13) 

C3Hs "IhT = 3 (r ^at + r DHPcat) < M) 

As the starting point for the kinetic analysis with the 
reactor model RPLUG from ASPEN PLUS, the most 
general rate equations for C0 2 shift reaction, FT reac- 
tion, and direct formation of hydrocarbons from C0 2 
(Table 2) were used. A first study with variable expo- 
nents for the partial pressures resulted in exponent 
values close to 1; therefore, the final modeling approach 
was done using a first-order law for all partial pressures. 
As a first approximation and in order to minimize the 
number of parameters to be determined, temperature 
effects on the adsorption constants were neglected in 
the temperature range investigated and the values 
estimated at 300 °C were used for all temperatures. 

Table 4 shows a summary of the results for the kinetic 
parameters which gave the best fit of the experimental 
data. For all rate equations, the best fit was obtained 
by considering inhibition terms according to the adsorp- 
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Table 4. Kinetic Parameter Values As Determined by 
Regression Analysis for CO2 Shift Reaction, FT Reaction, 
and Direct Hydrogenation of C0 2 (Table 2) [Catalyst 
100/13/10/10 Fe/Al 2 Os/Cu/K; T = 300, 330, and 360 °C, 
p = 1 MPa, (pHjpcoJin = 3, r mod = 0.042-21.4 g-s/cm 3 ] 

direct 

CO2 shift FT hydrogenation 

CLi HjO 65 ^3 ^ 

* i0 ,mol/(s-g-MPa) 1.51 x 10 7 2.58 x 10 8 39.6 
£ A ,kJ/mol 55 72 20 

tion of CO, H2O, and CO2. This is in agreement with 
Zimmerman and Bukur, 22 who, for the traditional FT 
synthesis with H2/CO, obtained the best fit of the 
experimental data with a similar catalyst also by 
considering the inhibiting effects of H2O and CO2 in the 
FT rate equation. The fact that in this study inhibition 
by CO2 has to be considered in all rate equations may 
be explained, besides by the high partial pressures of 
CO2 applied, in terms of the high potassium content of 
the catalyst. Addition of K increases the surface basicity 
of the catalyst 15 and therefore favors the adsorption of 
CO2, which has an acidic character. The temperature 
effects on the kinetic rate coefficients kj were assumed 
to follow the Arrhenius law (eq 15). 

kj = k jt0 exp(-E A /RT) (15) 

The apparent activation energy values obtained (55 
kJ/mol for the C0 2 shift reaction and of 72 kJ/mol for 
the FT reaction) are in the range of those reported in 
the literature (Tables 1 and 3). However, one has to keep 
in mind that in the literature the rate equations applied 
for the FT reaction as well as for the CO2 shift reaction 
were different from the ones used in this study. There- 
fore, a direct comparison is actually not possible, which 
might explain why in the case of the CO2 shift reaction 
the value for the apparent activation energy is on the 
lower end of the data reported. In the case of direct 
hydrocarbon formation from CO2, no published informa- 
tion is available about the kinetics. The resulting 
parameter values for the latter case are connected to 
the assumed structure of the rate equation (see Table 
2), which must be seen as a first approach. 

Taking into account only the reaction rate constants 
at 300 °C, the value for the FT reaction [70.4 mol/(s-g- 
MPa)] is more than 100 times higher than the one for 
the direct hydrocarbon formation [0.59 mol/(s-g*MPa)]. 
The reaction rate value rpH for reactor inlet conditions 
(300 °C, ph 2 = 0.75 MPa, and p C o 2 = 0.25 MPa) is 
equivalent to rpr for a CO partial pressure as low as 
0.0004pco2- These calculated values and the experimen- 
tal data obtained at the lowest residence time applied 
indicate that, for the catalyst used and for the reaction 
conditions applied, direct formation of hydrocarbons 
from CO2 plays only a very minor role compared to the 
direct formation via CO. For a more detailed and 
mechanistic clarification, especially for the extent to 
which the direct formation contributes to the overall 
hydrocarbon formation and for the determination of a 
more correct rate equation for hdh, further investigations 
at low residence time are required, which are currently 
performed in Korea at KRICT. 

The proposed reaction scheme together with the rate 
equations and kinetic parameters fits the experimental 
data with good accuracy. This is shown as an example 
in Figure 3 for CO2 conversion as well as yields of CO 




*C02.exp. . % 

Figure 7. Parity plot for all experimental and calculated CO2 
conversion values at 300, 330, and 360 °C. Catalyst: 100/13/10/ 
10 Fe/M 2 03/Cu/K. p = 1 MPa, (puJpcoX = 3, r mod « 0.042-21.4 
g*s/cm 3 . 

and hydrocarbons at a reaction temperature of 300 °C 
and in Figure 4 for the corresponding selectivity values. 
A parity plot for CO2 conversion (Figure 7) for all 
reaction temperatures applied reflects the good agree- 
ment of the experimental data with values calculated 
with the model. Except for the lowest residence time 
applied, the deviation of the calculated from the mea- 
sured values is in the range of less than ±10%. Also 
the selectivity data for CO and total hydrocarbons agree 
with calculated values within an error range of ±15%. 

6. Summary and Conclusion 

The effect of residence time and reaction temperature 
on CO2 hydrogenation was. studied in a fixed-bed mi- 
croreactor at 1 MPa [(pHjpco 2 )in = 31 using a promoted 
iron catalyst with the composition 100/13.3/10.7/8.9 Fe/ 
AI2O3/CU/K (mass ratios). The reaction temperature was 
varied between 300 and 400 °C and the residence time 
between 0.042 and 21.4 g-s/cm 3 . Product analysis was 
achieved by means of GC in connection with a TCD for 
the inorganic and a FID for the organic components. On 
the basis of the experimental results, a reaction scheme 
for CO2 hydrogenation was proposed. A kinetic model 
of the overall reaction, including equilibrium con- 
straints, was developed and applied using numerical 
features of ASPEN PLUS. 

The following conclusions can be drawn from the 
results. Chemical equilibrium calculations display that, 
for CO2 hydrogenation, in contrast to the traditional FT 
reaction with H2/CO, equilibrium limitations exist for 
CO2 conversion. Maximum CO2 conversion achievable 
from an equilibrium point of view exhibits a much lower 
value for the CO2 shift reaction alone than in a situation 
where also hydrocarbons are allowed to be present. 

The experimental results reflect a two-step reaction 
mechanism with CO2 reduction to CO as the primary 
step (equilibrium-limited with kinetic control) followed 
by subsequent CO hydrogenation to hydrocarbons (with 
kinetic control). This sequence of reactions allows for 
higher overall CO2 conversion than the CO2 shift 
reaction alone. 

Separation of both main reactions (CO2 shift and FT 
reaction) using two reactors in series would be disad- 
vantageous. For the endothermic CO2 shift reaction, the 
first reactor has to be operated at high temperatures to 
achieve a reasonable CO2 conversion. Even at about 500 
°C, less than 50% of C0 2 can be converted to CO from 
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an equilibrium point of view. For the second reactor (FT 
reaction), a much lower temperature is required, and 
therefore cooling of the reactants would be necessary. 

In the temperature range 300-360 °C, paraffins and 
a-olefins are formed as the main organic products of CO2 
hydrogenation with small amounts of oxygen-containing 
components. At 400 °C, significant amounts of carbon 
deposited on the catalyst are formed. Application of the 
ASF equation to the measured product selectivities 
leads to chain growth probability values a similar to 
those found for the traditional FT synthesis with H*j/ 
CO. This means that in CO2 hydrogenation practically 
the same organic products can be formed as if the FT 
reaction starts from CO. 

Experiments with varying residence times indicate 
that extrapolation of measured CO selectivities toward 
zero C0 2 conversion leads to values of hydrocarbon 
selectivity higher than zero. This agrees with a hypoth- 
esis to be found in the literature 19 that, besides the two- 
step reaction via CO, also a direct formation of hydro- 
carbons from CO2 can take place. However, presently 
this cannot be seen as evidence that in the intrinsic 
reaction mechanism a direct formation of hydrocarbons 
from CO2 on one active site can occur. 

Kinetic modeling of CO2 hydrogenation was carried 
out using the reaction scheme proposed as a conclusion 
from the experimental results. With a PFR model and 
kinetic parameters determined by nonlinear regression 
analysis, the experimental data could be fitted with an 
error of less than 15%. Values found for the apparent 
activation energies of C0 2 shift (55 kJ/mol) and FT 
reaction (72 kJ/mol) are in the range of those reported 
in the literature. Comparing the resulting reaction rates 
for hydrocarbon formation, it was found that compared 
to the formation via CO as intermediate direct formation 
from CO2 plays only a minor role. The results obtained 
for this reaction appear to be preliminary; additional 
work is required. Altogether, with the numerical meth- 
ods presented, it is possible to model the reaction 
kinetics of C0 2 hydrogenation on a K-promoted iron 
catalyst in the PFR used. The resulting kinetics should 
be regarded as a formal description for the catalytic 
system investigated in this study. 
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Nomenclature 

a«.H 2 o = inhibition coefficient for H 2 0 in reaction i 
av = average 

a = chain growth probability 

6;.co 2 = inhibition coefficient for C0 2 in reaction i 

BPR = back pressure regulator 

calc = calculated 

dp = particle diameter, //m 

d r = reactor diameter, cm 

Ar//°3oo*c = heat of reaction at standard pressure, kJ/mol 
A R G = Gibbs free energy, kJ/rnol 
exp = experimental 

E a = apparent activation energy, kJ/mol 

/min = function to be minimized (eq 10) 

Fi = axial molar flux of component t, mol/(cm 2 -s) 

<t> = variable in eq 10 



kj = reaction rate constant for reaction j, mol/(s*g*MPa) 
kj.6 = preexponential factor in Arrhenius law (eq 15), mol/ 
(s-g*MPa) 

K'cq = equilibrium constant of CO shift reaction (eq 5) 
tfeq = equilibrium constant of C0 2 shift reaction C-X/JC^) 
i T = reactor length, cm 
tficat = total mass of dried catalyst, g 
«Hc.i = molar flow of component 1, mol/min 
Ncj = number of carbon atoms in component i 
Ndata — number of values measured in each experiment (ea 
10) 

Nexp = number of experiments in data set (eq 10) 
v,^ = stoichiometric coefficient 
Pi = partial pressure of component t, MPa 
''shi = reaction rate of CO shift reaction, mol/(s-g) 
rsH2 = reaction rate of C0 2 shift reaction, mol/(s-g) 
r D H = reaction rate of direct hydrocarbon formation, mol/ 
(s-g) 

rpr = reaction rate of FT reaction, mol/(s-g) 

Peat = catalyst bed density, g/cm 3 

R = gas constant = 8.314, J/(mol-K) 

Si = selectivity of component i, C % 

o = standard deviation for measured values (eq 10) 

t = time, s 

Tmod = modified residence time (eq 7), g*s/cm 3 

T = temperature, K, °C 

y.o. = volatile organic products 

Vi r= volume flow at 0 °C and 1 bar, cm 3 /min 

x = reactor axial length coordinate, cm 

I, = extent of reaction y, mol 

Xco 2 = conversion of C0 2 , % 

Y, = yield of component t, C % 
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